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Structure of polymer electrolytes membranes, 
e.g. Naflon, inside fuel cell catalyst layers has 
significant impact on the electrochemical activ- 
ity and transport phenomena that determine cell 
performance. In those regions, Nafion can be 
found as an ultra-thin film, coating the cata- 
lyst and the catalyst support surfaces. The im- 
pact of the hydrophilic/hydrophobic character of 
these surfaces on the structural formation of the 
films and, in turn, on transport properties, has 
not been sufficiently explored yet. Here, we re- 
port about Molecular Dynamics simulations of 
hydrated Nafion thin-films in contact with un- 
structured supports, characterized by their global 
wetting properties only. We have investigated 
structure and transport in different regions of the 
film and found evidences of strongly heteroge- 
neous behaviour. We speculate about the impli- 
cations of our work on the understanding of real 
systems. 

Introduction.- Polymer electrolyte fuel cells [H-dl are 
complex devices constituted by soft and hard materials, 
which are characterized by heterogeneous nano-structure 
and transport properties. In particular, the membrane- 
electrode assembly has a strong impact on the fuel cell 
performance. It is formed, together with other compo- 
nents of limited interest here, by the polymer electrolyte 
membrane and the catalyst layer. The membrane must 
be an insulator against electrons transport, promote pro- 
ton conductivity also in low hydration and high tem- 
perature conditions, and be mechanically stable. The 
most successful material used in these days is Nafion 
which, upon hydration, phase-separates in hydrophobic 
and hydrophilic domains. The former provide mechanical 
strength, while the latter result in regions available for 
transport of water and protons. While the macroscopic 
structure of Nafion has been probed in depth, some de- 
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tails of the nano-structure still remain controversial. 

The catalyst layer is the interface region, whose com- 
plexity arises from the fabrication process. Here, carbon 
supported platinum nanoparticles mix with the ionomer 
and self-organize into extremely heterogeneous struc- 
tures. Despite significant efforts in the understanding 
of the bulk Nafion morphology, only a few studies have 
been devoted to the structure of the hydrated ionomer 
at the interfaces with platinum and carbon. Indeed, the 
ionomer is expected to self-organize in different forms, 
depending on substrate properties such as chemical com- 
position, geometry and, ultimately, wetting behaviour. 
The features of the film have significant impact on many 
processes including, amoiig others, the electrochemical 
double layer structure [1, |7| with the associated charge 
distribution, as well as proton conductivity Q. 

Computer simulations have been used for a better un- 
derstanding of the catalyst layer structure. Studies rele- 
vant in the present context include coarse-grained Molec- 
ular Dynamics simulations of the self-organized struc- 
ture of the catalyst layer [13] ; atomistic modelling of 
a) ionomer and water adsorption on graphitized carbon 
sheets ll'] and b) interfaces where the ionomer is in con- 
tact with the vapour phase and with the catalyst or the 
catalyst-support surfaces [l3] ; simulation of the effect of 
hydration of Nafion in the presence of platinum nanopar- 
ticles All these works are based on very detailed 
modelling of the solid phase in contact with the hydrated 
ionomer, both at the level of nano-structure and hetero- 
geneity of interactions. 

Here we have chosen a different approach, which fo- 
cuses very little on chemical details but aims at grasp- 
ing the main general physical features. This strategy 
has been successfully applied in studies of molecular 
liquids, like pure water in contact with infinite struc- 
tured walls [ij], or even simple atomic supercooled liq- 
uids [HI, [l^ . We have considered a mean-field- like inter- 
action ionomer/support, that allows us to precisely con- 
trol the hydrophilic character of the substrate by using 
a unique tunable control parameter. By using Molecu- 
lar Dynamics simulations, we have explored self-assembly 
and thin- film formation of a detailed model for the Nafion 
ionomer, when placed in contact with unstructured in- 
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FIG. 1. Snapshots of hydrated Nafion ultra-thin films, 
at T = 350 K and A — 22, for an interaction with 
the support of increasing hydrophilic character {e^ ~ 
0.25, 1.0, 2.0 kcal/mole, from top to bottom). We observe the 
formation of extended water pools (red) which are separated 
from the confining polymer matrix (grey) by the charged sul- 
fonic groups interface (green). For = 2.0 kcal/mole the 
ionomer is completely desorbed from the substrate. 



finite hydrophobic and hydrophilic surfaces (supports), 
under high hydration conditions. Also, we have focused 
on transport properties, and found evidences of strongly 
heterogeneous behaviour in different regions of the film. 
In what follows we report our main findings. Further 
results and a more comprehensive discussion will be in- 
cluded in forthcoming publications [U [3 ■ 

Model.- The Nafion ionomer is composed by a hy- 
drophobic poly-tetrafluoroethylene backbone with inter- 
calated perfluorinated side chains, which are terminated 
by a hydrophilic 5*03" group. An united-atom descrip- 
tion is used for representing C-P^.s groups in the polymer 
chain. In contrast, radical sulfonic acids, water molecules 
and hydronium ions are described with all-atoms res- 
olution. This hybrid modelling scheme is largely used 
to represent complex polymer systems. Intra-molecular 
bonded interactions include bonds, angles, and dihedrals 
terms. Non-bonded interactions are described in terms 
of Lennard-Jones potentials and long-range Coulombic 
interactions. Potential parameters are similar to those 
of Refs. [1^ [131 ■ Water molecules are described by the 
rigid SPC/E model [2l|, while the flexible model for the 
hydronium ion has been taken from Ref. ^2^. 

All system units interact with a smooth unstructured 



FIG. 2. Mass probability distributions as a function of the 
distance from the support, z, at the indicated values for e^j. 
We have considered a) water oxygens , b) polymer backbone 
units, and c) sulfur atoms. 

wall, placed at z = and parallel to the xy-plane, via a 
9-3 Lennard Jones potential. This only depends on the 
distance, z, of the imit from the support: 




-^), (1) 



where Zc = 15 A is a cut-off distance and 9 is the 
Heaviside function. The index a identifies complexes 
{H2O, H^Q'^, SO^) with significant dipolar coupling 
to the (hydrophilic) support (a = phyl), or units cor- 
responding to the hydrophobic sections of the polymer 
(a = phob) which, in contrast, interact very mildly. The 
energy well e^°^ — 0.25 kcal/mole, while eg,^^' = is 
our control parameter. The typical interaction length 
scale cr^j = 3. 2 A in all cases. We have calculated for 
Cw € [0.25, 2] kcal/mole the values of the contact angle 
of a droplet of water molecules interacting with the sup- 
port via Eq. ([T]), and demonstrated that this range en- 
compasses hydrophobic to strongly hydrophilic behaviour 
(lower to higher e^, respectively) ilS]. 

We have simulated systems formed by 20 polymers 
chains, containing 10 pendant side-chains each. Electro- 
neutrality was kept by adding 200 hydronium molecules. 
We have considered high hydration condition A = 22, 
where A is the number of water molecules per sulfonic 
acid group. Periodic boundary conditions were im- 
posed in the xy-plane containing the support, while open 
boundaries were kept in the z-direction. We have fixed 
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a temperature T = 350 K, and determined the box size 
(Lj. = Ly o± 80A) by matching a pressure P ~ atm. 
Initial system configurations were prepared by randomly 
placing ionomer, hydronium ions, and water molecules in 
the simulation box and letting all units self-organize in 
contact with the support. Less conventional computing 
details include calculation of long-range coulombic inter- 
action in slab geometry; thermostat coupled to the xy- 
components of velocities only; detection and deletion in 
subsequent data analysis of molecules evaporated during 
the molecular dynamics trajectory (a molecule is con- 
sidered as evaporated, and therefore disregarded, if at 
any time the z-coordinate exceeds 50 A, the typical thin 
film thickness). For production runs we have used the 
large-scale parallel Molecular Dynamics simulation tool 
LAMMPS dl. 

Results.- In Fig. [T] we show typical snapshots of the 
ionomer ultra-thin film for interactions with the sup- 
port of increasing hydrophilic character, from top to 
bottom. We observe the hydrophobic character of the 
ionomer /air interface and the formation of extended hy- 
drophilic pools, separated from the confining polymer 
matrix by an interface decorated by the sulfonate groups. 
Significant re-organization of the ionomer on increasing 
Cw is also visible. In order to quantify these changes, we 
plot in Fig. [5] the normalized mass probability distribu- 
tions for a) water, b) polymer backbone and c) sulfur 
atoms. For the most hydrophobic cases we have found a 
sandwich structure, with a wide water distribution cen- 
tred at z ~ 20 A and most part of the polymer in contact 
with the support and at the open interface with air. Dis- 
tribution of the sulfonate groups coherently shows two 
peaks at the ionomer/water interfaces. By increasing e^,, 
the ionomer abruptly desorbs from the support, induc- 
ing a significant re-organization of water domains. Water 
completely floods the region close to the support, where 
density layering is observed. A substantial decrease of 
hydration in regions far from the support is an expected 
consequence. 

From the above observations we can predict an im- 
pact of changes in morphology on local mass transport 
properties of water molecules in different regions within 
the film. These are, indeed, expected to depend on the 
interplay between confinement due to the hydrophobic 
matrix and direct interaction with the sulfonate groups. 
To quantify space-dependent transport, we restrict to 
the xy-plane, slice the film in overlapping slabs of width 
5z = A A equally spaced by lA, and consider the following 
generalized form of the mean-squared displacement: 

1 ^ 

(r2(i,z)) = -^(|r,(t)-r,(0)p5(z.(0)-z)). (2) 

Here () is the thermodynamic average and ri{t) is the 
xy projection of the three-dimensional position vector of 
the oxygen atom of water molecule i — 1, . . . ,N. Eq. ([2]) 
corresponds to the dynamics of molecules which, at time 
i = 0, were at a distance z away from the support. We 
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FIG. 3. a) Water molecules diffusion coefficient at the indi- 
cated values of e^,. Inset: Total diffusion coefficient, D, (cir- 
cles) integrated on the entire film. D is almost constant on 
the investigated range, but always higher than the value cal- 
culated in the membrane (dashed line). Main panel: D{z), as 
calculated from Eq. b) z— dependent structural relaxation 
time for water molecules, extracted from the self-intermediate 
scattering function. 



are now in the position to define a z— dependent diffu- 
sion coefficient, D{z), via the Einstein relation. This 
strategy has been demonstrated to properly characterize 
space-dependent translational diffusion in confined ge- 
ometries [T5l.[l6|. 

We first show in the inset of Fig. EJa) the total dif- 
fusion coefficient for water molecules, integrated in the 
entire film, compared to the analogous quantity calcu- 
lated in the membrane, far from any boundaries. (The 
diffusion coefficient of our pure water model in the same 
conditions is ~ 7.9 x 10~^ cm^/s.) Interestingly, we ob- 
serve a diffusion in the film augmented (of about 15%) 
compared to the membrane, at all values of e^,. Note 
that, although some kind of modulation with e^, is ob- 
served, we cannot conclude any continuous dependence of 
diffusion on the degree of hydrophilicity of the support. 
Next, we plot our data as a function of z in the main 
panel of Fig. [2Ua). For e^j^, < 1 kcal/mole, D{z) is quite 
uniform across the film, and everywhere close to the av- 
erage total value. This is consistent with a picture where, 
at high hydration, water is embedded in a crowded con- 
fining environment which is anyhow quite homogeneous 
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FIG. 4. z— dependent rotational relaxation time, extracted 
from an appropriate time-correlation function of the unit vec- 
tors normal to the water molecular plane. It is clear a strong 
correlation with the sulfonate groups mass distribution of 
Fig. He). 



within the film. In contrast, for > 1 kcal/mole dra- 
matic changes are visible, following the accumulation of 
the ionomer at intermediate and large distances, with 
the consequent formation of pure water layers at small 
z. Diffusion becomes strongly heterogeneous, steadily 
decreasing across the film with a rate which increases 
with e^. In the region close to the support, diffusion is 
enhanced compared to the integrated value, as also re- 
ported in the case of simple liquids confined by smooth 
boundaries [l^. At intermediate distances, D{z) is close 
to the integrated values, becoming strongly suppressed 
at higher distances (z > 25 A), where interactions with 
polymer backbones and sulfonate groups are very strong. 

We have checked the validity of these re- 
sults by extracting a space-dependent struc- 
tural relaxation time, t{z), from a general- 
ized one-particle intermediate scattering function, 
= l/7VE,(expiq.[r,(t)-r,(0)],5(z,(0) - z)), 
where q is the wave- vector. We have considered the 
smallest 2-dimensional q compatible with the simulation 
box, |q| = 2tt/Lx ~ O.OSA^^. Our findings are shown 
in Fig. [3l^b) and support our interpretation of diffusion 
data. We note that the apparent inversion of the 
two curves at largest e^'s can be ascribed to difficult 
sampling at large z, with no qualitative impact on our 
conclusions. 

Even more interesting is the effect of film nano- 



structuration on orientational dynamics. In Fig. |4] 
we show the rotational relaxation time, Tfl(z), ex- 
tracted from the correlation function C2{t,z) — 
1/^ (^2 (u, (t) • u, (0) )^ (z, (0) - z)} . u, is the vector nor- 
mal to the molecular plane of water molecule i, and P2{x) 
are the second order Legendre polynomials. Here the ef- 
fect of the interaction of water molecules dipoles with 
sulfonate groups is striking, with a non-monotonic mod- 
ulation of the rotational relaxation time, which strictly 
follows the 5*0;^ mass distributions of Fig. [^Jc). 

Discussion.- We have studied by Molecular Dynam- 
ics simulation the formation of Nafion thin-films at the 
interface with unstructured supports, characterized by 
their global wetting properties only. By tuning a single 
control parameter, e^,, we have been able to investigate 
in a unique framework a variety of environments pecu- 
liar of the catalyst layer, ranging from hydrophobic to 
hydrophilic (representative of carbon and platinum, re- 
spectively). We have found that an increase of the hy- 
drophilic character of the support has strong impact on 
the conformation of the film, which transforms from an 
irregular lamellar phase (where an extended water pool 
is sandwiched by ionomer sheets) to a phase-separated 
configuration, where water floods the interface with the 
support and polymers accumulate at the top. Intrigu- 
ingly, these mutations do not have strong impact on av- 
erage transport features of water molecules, as measured 
from the diffusion coefficient integrated on the entire film. 
In contrast, by implying a more involved analysis of dif- 
fusion at different distances from the support, we have 
discovered that dynamics at large is highly heteroge- 
neous across the film. We have shown that the diffusion 
coefficient is enhanced above the average value in regions 
close to the support, and strongly suppressed close to 
the ionomer /air interface, with a rate continuously con- 
trolled by e^. Moreover, very close correlation exists be- 
tween rotational dynamics and spatial sulfonate groups 
distribution. 

In conclusion, we intend to check if the observed 
features of the local environment embedding water 
molecules at high e^, , at distances from the support where 
layering effects are negligible, are similar to those present 
in the highly hydrated bulk membrane, far from any 
boundaries. In this case, our data would provide the first 
direct evidence of heterogeneous water transport in hy- 
drophilic channels as a function of the distance from the 
ionomer matrix boundaries. We are also convinced that 
our observations can be relevant in view of optimisation 
of real fuel cell catalyst layers. 
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